Color-charge separation in trapped SU(3) fermionic atoms 
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Cold fermionic atoms with three different hyperfine states with SU(3) symmetry confined in one- 
dimensional optical lattices show color-charge separation, generalizing the conventional spin charge 
separation for interacting SU(2) fermions in one dimension. Through time-dependent DMRG simu- 
lations, we explore the features of this phenomenon for a generalized SU(3) Hubbard Hamiltonian. 
In our numerical simulations of finite size systems, we observe different velocities of the charge and 
color degrees of freedom when a Gaussian wave packet or a charge (color) density response to a local 
perturbation is evolved. The differences between attractive and repulsive interactions are explored 
and we note that neither a small anisotropy of the interaction, breaking the SU(3) symmetry, nor 
the filling impedes the basic observation of these effects. 
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One of the most intriguing effects of strong correlations 
in low-dimensional systems is the separation of charge 
and spin degrees of freedom. As a generic phenomenon, 
a quantum particle carrying spin and charge converts 
in separate spin (spinon) and charge (holon) excitations 
with generally different velocities. On a microscopic bare 
level, there are rare examples where this process can be 
exactly studied in terms of explicit {n species) spinon and 
holon wave functions, as is the case for the Kuramoto- 
Yokoyama model [l|-|4j. However, on a low-energy ef- 
fective field theory level, spin-charge separation (SCS) 
explicitly manifests itself in all generic one-dimensional 
interacting systems belonging to the Luttinger liquid uni- 
versality class, where spinons and holons are described 
by independent collective density- type excitations [B| . In 
spite of several attempts, the direct experimental obser- 
vation of SCS has proved elusive So far, the best 
experimental evidence is provided by tunneling between 
quantum wires where interference effects relate to the ex- 
istence of excitations with different velocities ■ 

Since recently, one possible avenue for the observation 
and exploration of SCS are ultracold Bose or Fermi gases 
confined in optical lattices that have become an impor- 
tant instrument for investigating the physics of strong 
correlations. The great advantage of these systems is 
that the interaction parameters and dimensionality can 
be tuned with very high precision and control by means 
of an atomic Feshbach resonance or by changing the 
depth of the wells in an optical lattice [ll]]. One of 
the first achievements in this subject is the experimental 
observation of a superfluid to Mott insulator transition 
in a three-dimensional optical lattice with bosonic 87 Rb 
Interesting experimental results have also 



the observation of a Mott insulator in an optical lattice 
EH. Several theoretical works have addressed the_pos- 
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atoms 

been obtained for fermions, like the study of the molecule 
formation through the Feshbach resonance [l3|, the ex- 
perimental observation of fermionic superfluidity fbj] . and 



sibility of observing SCS in cold-fermionic gases 
and in cold-bosonic gases [22| with two internal degrees of 
freedom. Thanks to the special properties of cold atomic 
systems these theoretical proposals could address previ- 
ously unexplored features of SCS, like the problem of the 
breaking of SU(2) symmetry and SCS at high energies. 
In particular, higher spin fermions can be directly stud- 
ied with cold atoms in more than two hyperfine states. 
This kind of systems could give rise to new phases in opti- 
cal lattices due to the emergence of triplets and quartets 
(three or four fermion bound states), and other phenom- 
ena (23l - [29j . At least two alkali atoms 6 Li and 40 K seem 
to be possible candidates for the experimental realiza- 
tion of an SU(3) fermionic lattice with attractive interac- 
tions [2^]. In the case of 6 Li the scattering lengths for the 
three possible channels of the three lowest lying hyperfine 
levels (\F,m) = 1 1/2, 1/2), [1/2, —1/2), and |3/2, — 3/2» 
at large magnetic fields become similar for the three of 
them a s ~ — 2500ao Moreover, the realization of 

a stable and balanced three-component Fermi gas has 
been recently reported to potentially accomplish both an 
attractive and rep ulsive regime with approximate SU(3) 
symmetry [3l|,[32|]. The scattering lengths of the different 
channels for the three lowest hyperfine states of 40 K near 
the Feshbach resonance was also measured and the possi- 
bility of trapping them optically was demonstrated [33| . 

It is the purpose of this work to use time-dependent 
Density Matrix Renormalization Group (td-DMRG) (33 - 
liij l simulations to study the phenomenology of CCS in 
lattice systems with three different kinds of fermions, 
where the color notation is inherited from the quark de- 
scription of high energy SU(3) theories of quantum chro- 
modynamics. As one of the significative quantities ex- 
tractable from td-DMRG, the different color and charge 
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velocities are taken out from timc-dependent simulations, 
and compared to a low energy bosonization approach. 

The low energy physics of cold fermionic atoms with 
three different hyperfine states trapped in an optical lat- 
tice can be described by an SU(3) version of the Hubbard 
Hamiltonian, 



field 



H = -t 
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The sums a and /3 extend over the three colors rcd(r), 
green(g), and blue(b) corresponding to three different hy- 
perfine states. The operators fj a and f ia are the creation 
and destruction operators of an atom on site i with color 
a. We consider different values of the on-site interaction 
between the different color pairs U a p to be able to in- 
clude SU(3) symmetry breaking terms. The site label i 
goes from to L— 1, with L being the total number of lat- 
tice sites corresponding to the wells forming the optical 
lattice. For cold atoms, there is an additional harmonic 
confinement term that arises due to the Gaussian profile 
of the laser beams. If this trap potential is weak, we can 
assume to sit in the trap center where the confinement 
can be considered constant. L_arger potentials could be 
taken into account, e.g. 



see 
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In the subsequent cal- 
culations we will ignore the confinement. The hopping 
term can be controlled by varying the depth of the wells 
through the laser intensity. The optical lattice potential 
that each of the hyperfine states is affected by is very sim- 
ilar and the hopping rates can be considered to be equal 
for the three different colors. In the rest of the paper, all 
energies will be expressed in units oft — 1 . The interac- 
tion strength in each channel is proportional to the cor- 
responding s-wave scattering length U a p = Anh 2 a a p/m. 
The condition for the atoms to stay in the lowest band 
is that the s-wave scattering length must be smaller than 
the typical size of the wave function of an atom in one of 
the lattice wells [3!| . In turn, this is easily fulfilled in the 
experiments, so we will neglect the population of higher 
bands. Spin-flipping rates are usually smaller than the 
escape rate of atoms from the optical lattices and can be 
neglected as well. 

In the elementary linear bosonization approach, valid 
in the weak-coupling limit, the low-energy effective the- 
ory of the model can be expressed in terms of the collec- 
tive fluctuations of the densities of the three spin species 
plus charge density. Introducing the three bosonic fields 
4> a {x), the density operators for each color can be written 
as 



Pi, a ~ P+ -^=d x (f) a {x) - - sin 



2k F x + V47r^> Q (2;) , (1) 



where Uf is the Fermi wave- vector and x corresponds to 
the lattice site. We can express the bosonized Hamilto- 
nian in terms of the total density described by a bosonic 
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and the relative density described by two bosonic fields 



*3 = ^|(0r-0&) 
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The subindices were chosen to correspond to the SU(3) 
group Casimir operators J3 and Js- For a SU(3) sym- 
metric Hamiltonian, the model can be separated into two 
parts, charge and color, H = iJ c h + H co \, 



rr ^ch 



^(d x f> ch ) 2 + K(d x O ch ) 2 
is. 



(5) 



where u c h is the density velocity v c h = vf/K, K = (1 + 
2[//7rui?)~ 1 / 2 is the Luttinger liquid parameter, and Q c h 
is the conjugate field of the bosonic field $ c h, and 



H col = 

H=3,S 

u 
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■f — t COS V 27T$3 COS V 67T$ 8 — T — 7 COS V 87T$ 3 . 

(6) 

Similarly, the color velocity can be approximated as 
"col = vf \/l — U / (ttVf)- However, due to the non- linear 
cosine terms in (j6]), we do not expect a linear relation 
between distance and time for color excitations to hold 
for long times. If SU(3) symmetry is not strictly ob- 
served there appear mixing terms coupling density and 
color degrees of freedom. For example, if U rg = U r b = U\ 
and U g b = U2 the mixing term can be written as J? m i x = 
V2{Ui-U 2 )/3TT-d x <f> 3 d x <f>$ 0. A renormalization group 
analysis of this model can elucidate the low-energy prop- 
erties of our system [4(J Hl[ . The most important differ- 
ence with respect to the SU(2) case is that when U > 
but weak, umklapp processes present for commensurate 
fillings do not open a gap in the charge sector. We ex- 
pect a phase transition between the MI and the LL at 
a finite value of U. In fact, with Monte-Carlo calcu- 
lations, Assaraf et al. estimated the critical interaction 
scale at U CI = 2.2 [40| , while recent DMRG calculations 
suggest the critical point to be much closer to zero [42| . 
The cosine terms for the color Hamiltonian are irrelevant 
for repulsive, but relevant for attractive interaction and 
responsible of a gap opening in the color sector in the 
attractive case. 

Let us now elaborate on how the CCS appears in real- 
time simulations. We study the time evolution of the 
Hamiltonian (J) with the td-DMRG algorithm (15, [35|. 
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The number of states needed to keep sufficient accuracy 
during the time evolution was more than 7000. Such 
huge demands limited the system size we could simulate 
to L = 18 for periodic (PBC) and L = 48 for hard-wall 
boundary conditions (HWBC), which in total is settled in 
the range of present state of the art limits. For the small 
systems, the finite size effects are large. Comparable sim- 
ulations on SU(2) systems provided detailed knowledge 
about the finite size effects, so as that they are relatively 
independent on the interaction parameter. Comparison 
with exact results for U = of the charge and color ve- 
locities provide us with approximations to the Luttingcr 
liquid parameter. 

We show snapshots of the time evolution for differ- 
ent interaction strengths in Fig. [1] for the system with 
L = 18 and PBC. We calculate the ground state |\&o) of 
the SU(3) Hamiltonian with N = 27 fcrmions and then 
put an extra (green) fermion with an initial wave packet 
in our system that travels to the right with a finite mo- 
mentum centered around k = 0.6it 



\V+\t = 0)) 



Mkx) - 



fl\*o), (7) 



where the width of the wave packet is chosen to be 
a = 1.5. Choosing N ~ 27 leads to an average den- 
sity of (n) = 1.5 and an incommensurate filling of 
v = 0.5, the commensurate fillings being v = 1/3,2/3. 
Time is always measured in inverse units of the hopping 
rate. The panels in the figure show the particle density 
(n) = (ty +1 (t)\ J2 a fifa W) and the corresponding 
densities of the color quantum numbers j 3 and j 8 relative 
to the uniform ground state density at initial and a finite 
time. We observe a slight dispersion effect due to the 
finite width of the wave-packet. Moreover, the velocity 
of the excitation is not exactly the velocity of a plane 
wave with momentum k. However, as relative velocity 
quantities are concerned, without interaction the charge 
density and color velocity are exactly the same (upper 
left panel). In Fig. [TJ for the case U = 1, we start to see 
the separation of density and color degrees of freedom. 
The initial excitation separates and the color and charge 
density evolve with different velocities. The decay of the 
charge density excitation is more rapid for strong repul- 
sive interaction as it is seen in the case U = 5, clearly 
indicating the opening of the charge gap. In similar ways, 
a gap in the color sector should open for attractive in- 
teraction, and spoil the color density evolution in that 
regime. However, in our example with U = — 1, we can 
still observe rather clean and stable CCS, since, for our 
finite system, the arising gap in the color sector is too 
small to reasonably detect an enhanced color excitation 
decay. To extract color and charge velocities, we applied 
two different extraction methods. First, the velocity can 
be obtained dividing the number of sites traveled by the 
maximum of the density by the time. It is more accu- 
rate to fit a Gaussian on the density for every time step. 




FIG. 1: (color online) Snapshots of the time evolution of an 
additional fermion for filling v = 0.5 and U — (a), U = +1 
(b), U = — 1 (c) and U = +5 (d). The initial state is a 
Gaussian wave packet with average momentum k — 0.6n. In 
(b,c,d) the color density of quantum number j 8 (magenta) 
separates from the charge density (black), while j 3 (cyan) is 
constant by construction. Lines serve as guides to the eye. 
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FIG. 2: (color online) Snapshots of a perturbation at time 
T = (thin lines) and a finite time step (thick lines and data 
points). The dashed (solid) lines are fits to the color density 
of j 8 (charge density n) with the corresponding data points. 



There is a short transient time at the beginning, followed 
by a plateau of constant velocity until the packet hits the 
boundary. We extracted the velocity of the packet when 
it passes the median between the initial position and the 
boundary and took a Gaussian averaged value around 
this position. As a second measure, for a definite plateau 
of constant velocity, we used a L = 48 sites system with 
HWBCs and replaced the incident electron by an initial, 
small potential perturbation. The time development of 
this method differs mainly by the implicit stimulation of 
excited states around both Fermi points leading to sym- 
metric peaks running in both directions and increasing 
the transient time at the start. Fig. [2] shows snapshots 
of the initial and a finite time step for U = 1.5, where the 
perturbation was taken to be the derivative of a Gaus- 
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that high anisotropics totally destroy charge-color sepa- 
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FIG. 3: (color online) Velocity v(U) of an initial Gaussian 
perturbation. Results of our td-DMRG simulations for both 
charge density (+) and color density (x). The dashed lines 
represent bosonization results, scaled to match the DMRG 
velocity at U = 0. The solid line is a fit on a next-to linear 
dependency in U. 



sian and this function was fitted on the two propagating 
wave packets. We find that both methods provide sim- 
ilar results, while the latter proved to be a considerable 
improvement in precision and is put forward by us as a 
general preferable method to extract velocities from td- 
DMRG. 

Finally, the extracted velocities are shown as a func- 
tion of the uniform interaction in Fig. [3J On top of the 
data, we show the expected relation for v co \(U), v^iJI) 
from bosonization, scaled to match the free fermion ve- 
locity. This a priori renormalization covers all direct ef- 
fects that lead to deviations of the Fermi velocity due 
to finite width or the central momentum fco of the ex- 
citation and the finite size dispersion. Even for \U\ up 
to 1, the agreement between the bosonization estimate 
and the numerical data is found to be excellent (dashed 
lines in Fig. [3]). Beyond, the limits of the low energy 
expansion are apparently reached. However, upon ex- 
tending the velocity expansion to higher order in U, e.g. 
Vch = vpy/l + 2U/ {nvp) + all 2 , the fit already covers 
the complete computed range in U (solid line in Fig. [3]) 
and the parameter a = —0.017 has the same order and 
sign as the fit by Assaraf et al. [!(| who chose a different 
filling in their Monte-Carlo simulations. 

A possible candidate for realizing the trionic phase 
(U < 0) is 6 Li for which the magnetic field dependence 
of the three scattering lengths has been measured [3p| . 
The attractive U interaction for magnetic fields around 
1000 G can be estimated as U rg = Uo, U r b = 1.23C/ , 
and U g b = 1.06[/o, due to the ratios in the s-wave scat- 
tering rate in each channel j30l |. It is important to con- 
sider whether these anisotropics infringe on the valid- 
ity of the effects observed above. In bosonization the- 
ory, anisotropy introduces new terms that couple the 
charge and color degrees of freedom so we would expect 



ration 41] . However, small and experimentally accessi- 
ble anisotropics turn out to be not decisively important. 
We have checked that the separation effect of color and 
charge densities is still visible in the simulations without 
much change. Even for the special cases of commen- 
surate fillings where a stronger sensitivity on the color 
anisotropy may have been suspected, the qualitative be- 
havior persists. 

In summary, we have performed td-DMRG simula- 
tions of cold fermionic atoms with three hyperfine states 
trapped in an optical lattice. Our simulations allowed us 
to observe the color-charge separation in SU(3) fermionic 
systems in a generic non-commensurate case. We have 
obtained the charge and color velocities as a function 
of the interaction from the real-time simulations both 
for the attractive and repulsive case. Once we take 
into account finite size effects by renormalizing the non- 
interacting velocities, our results at weak coupling are in 
good agreement with bosonization calculations. 
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